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Using neutron scattering, we investigate the static and dynamic structure functions S(Q, cu) of 3 
distinct solid phases (using powder average techniques), and characterize the corresponding modes 
of energy excitation, focusing primarily on the dynamics of spin interaction. With the S(Q,uj) 
maps, we extract the temperature dependent production rate and upscattering cross section, that 
are relevant for understanding experimental data on ultracold neutron production in solid oxygen 
at the saturated vapor pressure. 
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I. INTRODUCTION 

Because of its non-zero electronic spin (5=1), oxygen 
displays a uniquely rich phase diagram when compared 
to other cryocrystals made of diatomic molecules, such as 
N2, H2 and CO. At the saturated vapor pressure, there 
exist 3 distinct solid phases: 7 (44 ~ 55 K, cubic lat- 
tice, paramagnetic), /3 (24 ~ 44 K, rhombohedral lattice, 
short range anti-ferromagnetic(AF)), and a phase ( < 
24 K, monoclinic lattice, 2-D AF)i. When subject to 
high pressures, oxygen loses it magnetism and this Mott 
insulator becomes metallic. Upon further pressurization, 
the solid becomes superconducting at temperatures be- 
low 0.2 K. While there exists great theoretical interest 
with the complexity of these phase transitions of this 
seemingly simple system, in the paper, we will focus on 
gaining understanding of the spin dynamics in the low 
temperature phases at low pressures, and illustrate their 
potential applications of ultracold neutron (UCN) pro- 
duction. 

The magnetic excitations in solid oxygen (S-O2) bring 
about an opportunity to enhance the superthermal UCN 
production through inelastic magnetic scattering, be- 
yond the typical phonon couplings already applied in 
supcrfluid helium and solid deuterium with a great suc- 
cess. In a superthermal UCN source, the dominant up- 
scattering loss is suppressed by reducing the number of 
thermally populated energy quanta (phonon or magnon) 
through cooling of the source material, whereas the pro- 
duction is relatively temperature-independent as long as 
the inelastic scattering channels remain open and unper- 
turbed. The phase space of the neutrons can be com- 
pressed through energy dissipation by exciting the collec- 
tive modes in the solid. In s-02, the magnetic scattering 
length is comparable to that of nuclear scattering, and its 
spin degree of freedom allows additional inelastic neutron 
scattering through magnetic coupling that has never been 
applied to neutron moderations. In this paper, we report 
results on modern inelastic neutron scattering (INS) in 



attempt to understand the physics of magnetic excita- 
tions in S-O2 relevant to our UCN production experiment. 

At low temperatures, the spin ordered state of S-O2 is 
a Neel state with a AF ordering in the basal plane. As 
the temperature rises, the spin-spin correlation remains 
strong enough to evade thermal fluctuations, resulting in 
a short-range spin order in the j3 phase, and the spin cor- 
relation persists in the 7 and liquid phases up to relatively 
high temperatures. The spin system does not follow the 
simple order-disorder transition, but instead an interme- 
diate phase exists between the spin-ordered a phase and 
spin-disordered 7 phase. This intermediate phase is best 
categorized as a geometrically frustrated spin system, a 
topic of which much work are actively ongoing to for- 
mulate viable wave-functions to calculate energy excita- 
tions^. Traditionally, the spin order in the j3 phase was 
described by a helimagnet or a three-sublattice Yaffct- 
Kittel structure^. However, due to its low spin number, 
the system should really be treated quantum mechani- 
cally with discrete spin quantizations. The classical spin 
models might not produce the correct excitation spec- 
trum^, especially in the strongly coupled, geometrically 
frustrated /3 phase. At low temperatures under zero pres- 
sure, the frustration is relieved by a structural distortion 
as the hexagonal lattice structure is slightly sheared into 
a monoclinic structure with reduced symmetry. 

The monoclinic unit cells of the a and j3 phases can be 
transformed into each other by a small inhomogeneous 
deformation of the closely packed basal planes, driven 
by the spin ordering. The c± lattice constant increases 
when approaching the phase transition from high temper- 
ature, leading to the dominant 2D spin alignment, with 
AF exchange interaction and an easy-axis anisotropy. 
The two-sublattice AF monoclinic lattice transforms dis- 
continuously into the three-sublattice AF hexagonal lat- 
tice and undergoes a second order phase transition into 
the magnetically disordered structure (albeit with short 
range spin correlations). It was suggested that magneto- 
elastic coupling changes the order of the transition from 
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the second to the first one with a sizable latent heatl. 
In addition, it was not yet understood why the a — j3 
transition happens at a temperature as low as 24 K, de- 
spite of its strong exchange field. Inter-plane exchange, 
spin-orbital, and magneto-dipole interactions all have the 
opposite effect to cut down fluctuations at the small-/s. 
This can probably be resolved if the quantum fluctua- 
tions in the quasi-two-dimensional Hcisenberg AFM are 
taken into account. 

Concerning the UCN production, we are most inter- 
ested in the energy excitation modes that incident cold 
neutrons (CN) can efficiently couple to. While the exci- 
tation in CI-O2 can be described by magnons using spin 
wave theory, the situation in j3-0— 2 still needs to be 
clarified. The first modern many-body calculation of ex- 
citations in S-O2 revealed interesting behaviors^. In the 
/? phase, the ground state has a two-fold degeneracy with 
a low lying first excited state at 0.8 meV. Because this 
excited state corresponds to the center of the Brillouin 
zone, as opposed to the zone boundary as with typical 
magnon excitations, the authors speculated the existence 
of bound pairs of spin waves in (3 phase. This hints to the 
formation of a spin singlet solid, a necessary condition to 
attain the spin liquid state. Another Monte-Carlo sim- 
ulation which allows a dcformable cell^ also showed an 
abrupt increase in spin dynamics in the j3 phase beyond 
that of the stable a phase. Comparing the excitation 
spectrum in S-O2 to other well studied frustrated spin 
systems, such as a-NaMn02 (with S = 2)2 and CS2CUCI4 
(S = 1/2 we find many similarities in the INS data. 



qucntly solidified and cooled under the saturated vapor 
pressure. A S-O2 sample in the polycrystalline form was 
attained under typical cooling procedures, as the ther- 
mal contraction, in particular the abrupt 5% change of 
the molar volume through the 7-/3 transition, would most 
likely crack the bulk solid into many micro-crystals. The 
cryostat was allowed to stabilize at each set temperature 
for 1 hour before beginning a measurement. Data were 
collected at 2 K temperature step throughout the (i and 
a phase with smaller temperature steps around the a-(3 
transition. In the final step, the cell was warmed up to 
totally evacuate oxygen and the background data were 
collected with the empty aluminum cell at 24 K. Cali- 
bration data on detector efficiencies were later collected 
using a vanadium foil (21.6 g in the beam). 

The data were reduced using the DCS MSLICE routine 
in the DAVE analysis suite developed by NCNR. Each 
data file was processed with proper background subtrac- 
tion and detector efficiency normalization. We used the 
powder average technique, in which signals from any de- 
tectors with the same 26 diffraction angle were averaged 
and the final results were presented as a function of \Q\ 
for each TOF bin (or w bin). The processed data were 
re-binned and exported as tables of S(Q,lu) for further 
analysis. 



III. FORM FACTORS FROM ELASTIC 
SCATTERING 



II. EXPERIMENTS 

We measured the total cross section of cold neutrons 
using a standard transmission measurement, as a part of 
the UCN production experiment at Lujan Center neu- 
tron flight path 12 at the Los Alamos Neutron Science 
Center (LANSCE) . Details of the experimental setup are 
described in — . The INS experiment was carried out 
using the Disk Chopper Spectrometer (DCS)^2 at the 
National Institute of Standards and Technology (NIST) 
Center for Neutron Research (NCNR), where neutrons 
are produced from a 20 MW fission reactor. Most of 
the data were taken using neutrons at 2.3 A (15.4 meV). 
With this wavelength, the energy resolution (FWHM) of 
the instrument is 1 meV at elastic scattering, and im- 
proves as the neutron time-of-flight (TOF) increases to a 
resolution of ~ 500 ficV at energy transfer u> = 15 meV. 
The target cell was mounted inside an "orange" helium- 
flow cryostat. Liquid oxygen was first condensed from 
gaseous O2 (at 99.999%) in an annular cell made of alu- 
minum, with an outer diameter of 12.5 mm and an in- 
ner diameter of 10.5 mm, leaving an annular thickness 
of 2 mm. With the height of 110 mm, the cell has a 
total volume of 6.1 cm 3 . The beam was collimatcd to 
90 mm of total height, and illuminated the entire width 
of the target cell. Inside the cell, the liquid was subse- 



In order to get a correct handle of the distinct form 
factors in magnetic and nuclear scattering, required later 
to calculate the inelastic cross sections, we start by ana- 
lyzing the data on elastic scattering. The total cross sec- 
tions for a- and P-O2 are plotted in Fig.Q]and the elastic 
diffraction lines are plotted in Fig. [2] Note that the to- 
tal cross section is predominantly elastic over the energy 
spectrum of the incident cold neutrons. The elastic cross 
section can be separated into coherent and incoherent 
contributions: 

< h = ^^yY,d hkl \F£$fe- 2W , (1) 

71 hkl 

a™ c = 4n\F mc \ 2 e~ 2W , (2) 

where dhki is the inter-planar spacing between the [hkl] 
planes in the lattice space, i.e., 

2ir 

dhkl = —=r, where G = hg x + kg 2 + lg 3 . (3) 
\G\ 

Here G is the reciprocal lattice vector with integer num- 
bers of [h,l,k]. The form factor F within the unit cell 
differs depending on whether the coherence of the scat- 



3 



Total Cross-section of Solid a-0 2 (monoclinic) 

18 . : : — 

nuclear, No FFs 
16 - - nuclear, isotropic FF 
nuclear, mole, axis FF 

^ nuclear, mole, axis FF + i//3 

.—. magnetic, AF, isotropic FF 

E • data (a-OJ 

oj 

5.10 magnetic AF, mag. FF 



[111] 

y 




(a) 



1£ 



20 
18 
16 
14 

CO 

c 

ffl 12 

§ 10 
o 

1> 8 

CO 

I 

co 6 

4 

2 



-2 
10 



10" 

Energy (meV) 



10 



Total Cross-section of Solid |3-0 2 (Rhombohedral) 



o 



nuclear, no FFs 
nuclear, isotropic FF 
nuclear, mole, axis FF, if/3 
data (P-OJ 



t t t t ▼ 



10 

Energy (meV) 



(b) 



FIG. 1: Total cross-section of a and /J-O2 
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In a non-Bravais crystal, the number of atoms in the unit 
cell, Nbasis can be more than 1. 

As shown in Fig. [T] form factors significantly modify 
the cross section in the following ways: For nuclear scat- 
tering, the intra-molecular interference between the two 
atomic nuclei would introduce an additional molecular 
form factor, 



pmol. 



JQ-r 



(6) 



where f is the vector representing the intra-molecular 
separation of the two atomic nuclei. To facilitate the 
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FIG. 2: Neutron diffraction lines of S-O2. The dashed lines 
are the 3 major background peaks from Al container cell are 
at Q=2.69, 3.11, and 4.37/ A. 



AF spin ordering, the O2 molecular axis is aligned to be 
perpendicular to the basal plane. Hence, the vector r is 
represented as (0,0,1), with I = 1.21 A. Note that it would 
not be appropriate to apply the typical spherical Bessel 
function used in diatomic molecules having free rotations, 
such as H2 and D^. Librational motions would lead 
to an additional Dcbyc- Waller factor. The inclusion of 
pmol. rec i UC es the elastic nuclear cross sections by about 
50%. 

For magnetic scattering, Stephens^ and Dunstetter— 
suggested the magnetic form factor to be used for molec- 
ular O2. In the a phase, we need to introduce a sign 
change in the unit-cell form factor in Eq. 2] to reflect on 
the AF spin configuration. In a conventional unit cell 
with two O2 molecules of opposite spin orientations, the 
lattice form factor becomes 



^coh 



b e 



1QR1 



b e 



iQB,2 



(7) 



where i?i = (0,0,0) and i? 2 =(a/2,b/2,0) arc the vec- 
tor of the two basis in a monoclinic unit cell. The 
lattice constants are a = 5.375 A, b = 3.425 A, and 
c = 4.242 A for mon oclinic a-02, and reduce to a = 
^3.754 2 + (3.425/2) 2 A for rhombohedral /3-0 2 . As the 
AF long-range spin order is present only in the a phase, 
the change of sign does not apply to the (3 phase. Dou- 
bling the periodicity of the spin lattice in a-02 mani- 
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fests as the additional Bragg edges at long neutron wave- 
lengths in the total cross section (at the E<1.5 meV in 

FigHk) and magnetic diffraction peaks below 1.7 A 
(Fig. [5k,). The lowest magnetic Bragg cutoff at 0.8 meV 
results from the constructive interference between the 
[100] planes, which is the perpendicular to the direction 
of spin order. The second lowest magnetic Bragg edge 
corresponds to scattering along the [101] plane. Due to 
its small intensity and close proximity to the first nuclear 
Bragg edge, this second magnetic Bragg edge is harder 
to resolve in the total cross section, but can be clearly 
identified in the powder diffraction data. The two low- 
est Q lines in the a phase (marked by the asterisks in 
Fig. [2k, first reported by Collins^ as magnetic in origin) 
disappear in the /3 phase, providing the evidence of long- 
range spin ordering in a-C>2- Also note the similarity of 
the total cross section and the rest of the diffraction lines 
between a and P-O2, indicating an adiabatic transfor- 
mation of the nuclear lattice structures with some broken 
symmetries in the a phase. We will take advantage of the 
similar lattice structures in a and P-O2 to separate the 
magnetic excitations from the phonon excitations. The 
7-O2 displays the symmetry of a cubic lattice (Fig. [3b), 
with quite distinct diffraction patterns. 



IV. DYNAMICS FROM INELASTIC 
SCATTERING 




The neutron scattering intensity is characterized by 
the double differential cross section, expressed by 



FIG. 3: The dynamic structure function S(Q,u>) of S-O2 in 
the (a) a phase at 4 K, (b) /3 phase at 32 K, and (c) 7 phase 
at 47 K, measured by INS. (d) Constant Q-slices of S(Q,u) 
in a-Oz (e) Q-dependent amplitudes of integrated S(Q, u>) in 
/8-Oa. 



For polycrystalline samples, the dynamic structure func- 
tion S{Q,u>) is measured as a powder average. In the 
limit that the orientations of micro-crystals are fully ran- 
domized inside the bulk sample, the angular dependence 
of the momentum transfer Q is averaged over 4-7T to be- 
come \Q\ = Q. During each scattering, the incident neu- 
tron experiences a momentum transfer of Q = k — k', 
and an energy transfer of w = Ei — Ef. Processes with 
positive uj involve neutron downscattering with the neu- 
tron losing energy into creations of one or more quasi- 
particles in the S-O2; processes with negative uj involves 
neutron upscattcring with the neutron gaining an energy 
gain by annihilating one or more quasi-particlcs, absorb- 
ing all their energies. Two dimensional maps of S(Q,uj) 
at different temperatures were produced with on aver- 
age 2 hours of neutron exposure. Samples of the S(Q, uj) 
maps are shown in Fig. [3] with color indicating the scat- 
tering intensity normalized by the phase space factor, 
k/k'. The INS measurements unfold the very distinct 
dynamics among the 3 solid phases. 



A. Q-O2 

As we have discussed in sec. IIII1 the long-range AF 
order is responsible for the two magnetic Bragg lines at 
Q = 1.32 A" 1 (Ml) and 1.57 A _1 (M2). To get a quanti- 
tative measure of the intensities of these elastic magnetic 
lines, we integrate the S(Q,uj) from uj = [—1,1] meV, 
and normalize to that of the first nuclear Bragg peak at 
Q = 1.7 A . The temperature dependence of the inten- 
sity of Ml can be compared with previous optical mea- 
surement^ (Fig. [U inset), and shows good agreements 
with the prediction of 2D Ising Model. The decreasing 
intensity with increasing temperature is consistent with 
the thermally induced reduction of the spin-spin correla- 
tion predicted in the Monte-Carlo simulation^. 

Following the work by Stephens et alJ±, we describe 
the magnetic excitations in a-02 with a Heiscnbcrg spin 
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magnetic clusters in a-02. Inset: Intensity of the magnetic 
elastic lines in a-02- 



Hamiltonian: 



H = -2 



E 

<ij> 



JijSi ' Sj 



-D^2s 2 xi -D'^2(S 2 yi -S 2 zi ), (9) 



based on the renormalized spin wave theory previously 
developed for NiCl?-^, which is a simple AF system with 
a rhombohedral structure. The exchange interactions in 
solid a-02 include the 4 nearest neighbors, Jnn, and 2 
next nearest neighbors, Jnnn, all within the basal plane. 
The inter planar coupling is negligible due to the large 
planar separation. The x axis is the easy magnetization 
direction (along the crystal b axis) , along which the spins 
align co-linearly. The y axis points along the crystal a 
axis, and perpendicular to the basal plane is the z axis, 
along which aligns the molecular axis. The Hamiltonian 
contains the single particle term with a self-energy D, 
introduced to characterize the spin anisotropy. The Ising- 
like behavior, as revealed by the temperature dependent 
intensity of the magnetic Bragg peak, strongly suggests 
that that the anisotropies shall satisfy D > and D' ~ 0. 
If D = D' > 0, we would expect the system to have an 
easy plane magnetization. 

Using the renormalized spin wave theory, the magnon 
energy spectrum of this two-dimensional S = 1 system is 
solved-----: 

E 2 a = (2J NNN {0) - 2J NNN (q) - 2J NN {0) + Df 

-[2J NN {q) + {-l) a D']\ (10) 

with the index a — as the acoustic mode and a = 1 
the optic mode. The coupling Jnn(q) = J2i JNN^' Ti 
is summed over the nearest neighbors on the opposite 
sublattice, and Jnnn is the next nearest neighbors on 



the same sublattice. The best value of the coupling con- 
stants reported by Stephens--' are: Jnn = —2.44 meV, 
Jnnn = -1-22 mcV, J± = 0, D = 0.134 meV, and 
D' = 0.120 meV, which suggests an easy-plane magne- 
tization (a XY- model). As we will show next, the cou- 
pling constants should be modified, in order to bring the 
model in better agreements with the new high resolution 
INS data and to conform with the Ising model suggested 
by the elastic scattering data. Fig. [5^ illustrates the en- 
ergy dispersion relations calculated using the updated 
coupling constants along major symmetry axes. Not sur- 
prisingly, there exists a finite energy gap for magnon ex- 
citations even at Q = 0, as a result of the anisotropy. As 
we allow the scattering to be confined in the basal plane, 
the S(Q, u>) starts to cluster around regions with modes 
crossing (Fig. [SJd) along the [100] direction. 

As neutron scattering probes fluctuations perpendic- 
ular to the scattering vector Q, the magnetic dynamic 
structure function of a AF system can be written as 

S ma °(Q,Lu) = (1 + Ql)[u q + (-l) m v q } 2 6(E a (q) - hio) 

xS(Q-q-G). (11) 

Here the scattering amplitude is enhanced when the mo- 
mentum transfer points along the spin aligned axis, i.e. 
the b axis. The delta functions enforce the energy and 
momentum conservation in this coherent scattering pro- 
cess. The Holstcin-Primakoff transformation coefficient 
used in Stephens' paper—, when compared against the 
experimental data of 5(Q, ui), is found to carry an incor- 
rect sign, and should be corrected as the following: 



*, + (-!)' 
,(2Jo 



2S- 



2Jo 



2J q + D) + (~l) m [2J' + {-l) a D'\ 



EM 



(12) 



It was this sign error together with the low values of 
anisotropy energies that led to predictions of magnon 
modes with energies lower than the result of the new 
INS measurements. This led us to over-predict the UCN 
yields when we first proposed the idea of magnetic UCN 
production using 

With the Q-dcpendcnt energy dispersion relation in 
Eq. [TU1 we calculate the dynamics structure function for 
magnon excitations in polycrystalline a-02, following the 
algorithm we have developed to calculate coherent in- 
elastic scattering cross section in solid deuterium-^. The 
magnetic form factors on the molecular and lattice level 
have already been validated using the elastic cross sec- 
tion as discussed in Sec. lIIII With the coupling constants 
suggested by Stephens et.al^l, the calculated S ma9 (Q, u>) 
is shown as an intensity plot in Fig. [SJ;. Because of the 
molecular spin form factor, the scattering amplitude be- 
comes negligible for Q > 2A , above which the phonon 
scattering dominates. The small anisotropy energies lead 
to a soft spectrum (clustered around the center of the sec- 
ond Brillouin zone) extending down to energies as low as 
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FIG. 5: Calculated S(Q,iS) of coherent magnetic scattering 
in a-02 using the magnon energy dispersion relation derived 
from the normalized spin wave theory, (a) along symmetry 
axes (b) inside the basal plane (c) polycrystalline average with 
small anisotropy energies (d) polycrystalline average with a 
large anisotropy energy. 



izes the inter molecular spin-orbit and spin-spin coupling. 
The single spin anisotropy term is A=3.96 cm -1 . When 
considering the molecular librations, the spin-libration 
interaction can only soften the anisotropy by 10%i. 

Varying the value of Jnn, on the other hand, shifts 
the energy of the 19 mcV magnons (corresponding to 
excitations at the Brillouin zone boundary), which is be- 
yond the range of energy coverage of the measurement. 
The difference Jnn — Jnnn determines the placement of 
the dominant magnon at 10 meV at the zone boundary. 
Results of our numerical studies agree with Stephens' 
suggestion of Jnn — Jnnn = —1-22 meV. 

We calculate the density of states of the magnetic ex- 
citations by inverting the dynamic structure function as 
detailed in later discussions (sec.[V]). As shown in Fig- El 
increasing the easy-axis anisotropy energy significantly 
shifts the population of low energy magnon modes to- 
wards a higher energy gap, whereas the 10 meV mode 
shifts only by 5%. The density of states derived from 
experimental data (limited to magnetic excitations at 

Q < 1.7 A ) is plotted for comparison. Because of the 
limited range of coverage in the energy transfer at low Q, 
the data cuts off around 12 meV. Except for the sharp 
peak around 10 meV (which is washed out by the fi- 
nite energy resolution no better than 0.5 meV), the mea- 
sured density of state of magnetic excitations is in general 
agreement with the results of the 2D Ising AF spin- wave 
model, with a large anisotropy energy. 

For Q = 1.2 ~ 1.4 A , the scattering is predominantly 
magnetic. We fit S(Q,ui) over the constant Q slice as a 
sum of a Lorentzian (centered around 10 meV) peak and 
a phenomcnological 3-polc response function Eq. [15] (as 
shown in Fig. As the solid temperature increases, 

both the high energy cluster (E^IO meV) and the low en- 
ergy cluster (E^5 meV) shift downwards in energy (see 
Fig. [H in accord with the thermally induced reduction in 
exchange couplings, as previously shown with the inten- 
sity of the elastic magnetic lines (Fig. 2] inset). 

The Neel temperature T N = lS{S + 1)J(0) = 
|(4Jjvjv — 2 Jnnn), calculated within the mean field ap- 
proximation, is 35 K. This is about 1.5 times higher than 
the a — j3 transition temperature of 24 K. Inclusions of 
the molecular field, however, could not explain the dis- 
crepancy. Perhaps the quantum fluctuations prevent the 
spin from ordering until the lower temperature is reached. 
The exchange field can be estimated using 



~ 0.8 meV. If D = 0, the excitation would become gap- 
less at Q = 1.3A . If D = D' > 0, the excitation also 
becomes gapless. It requires the adjustment of the easy 
axis anisotropy energy up to D = 0.6 meV, to bring the 
low energy mode higher than 5 meV to better agree with 
the experimental data (Fig. [3^,). This again supports the 
conclusion that the nature of the magnetism in a-02 is 
indeed 2D Ising. The increased anisotropy energy also 
agrees with the free gas value of 5.712 K, that character- 



2zJ N , 

Oexch. — — — -L 



(13) 



g^BS 

is quite large. With the measured spin flop field of 7 T— 



^sf ^J^anis.J^&xch.i 



we arrive at the anisotropy field of about 0.11 T. 
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P-O2 



As revealed with the INS measurements (Fig. [3]), P-O2 
has a magnetic excitation that differs significantly from 
the magnons in the a phase. First, the intensity of the 
magnetic Bragg line at Q = 1.3 A drops by orders 
of magnitude when compared to that in the a phase. 
However, there appears a prominent magnetic excitation 
centering broadly around Q = 1.3 A . The excitation 
is dispcrsion-less in energy and extends over lu = ~ 
8 meV. When compared to the spin diffusion dynamics in 

7-O2 at low Q < 1.5 A — , this dispersion-less magnetic 
excitation in (3-02 looks intrinsically different as it docs 
not extend continuously to low Qs as in 7-O2. There 
are many attempts to derive the spin configuration and 
its short-range order by studying the detailed structure 
of the broad elastic peak, based on different variations 
of the three-sublattice Yaffet-Kittel structur o 22 ' 23 . We 
would like to point out that the spin-pairing to form a 
spin singlet could also account for the vanishingly small 
clastic scattering amplitude. However, high quality data 
using single crystals are required to definitely support the 
spinon formation in /3-02- 

Using the dissipation-fluctuation theorem, the ampli- 
tude of inelastic scattering in a magnetic system can be 
described as the imaginary part of the magnetic suscep- 
tibility x{l) = x'(l) + *x"(<z)- To characterize the en- 
ergy dependence of this scattering, we fit the scattering 

amplitude (summed over Q = 1.2 ~ 1.4A ) using the 
phenomcnological 3-polc response function, 



, X'(g) 
2 

r(?) 



ujx 



r(«) 



(15) 



(c- Wo ( £2 )) 2 +r 2 ( (7 ) (LO + o Jo ( q )r + T^q) / 
where x'(o) is the q-dependent susceptibility, T(q) is the 
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FIG. 7: Inelastic magnetic scattering in P-O2 fitted with the 
phenomenological 3-pole response function (a), (b) The tem- 
perature dependence of the fitted parameters, (c) The tem- 
perature scaling of x'il)- 



inter-site fluctuation rate, and coo(q) is a phenomenologi- 
cal energy pole. Examples of the fit and the temperature 
dependence of the fit parameters are shown in Fig. [7] 

Most of the early literatures used the "coupled para- 
magnet" model to explain the large Lorentzian type of 
dynamic response and the g-dependence. However, the 
rhombohcdral structure in P-O2 is equipped with all nec- 
essary conditions of a geometrically frustrated AF mag- 
net. We attribute the observed dispersion-less energy 
excitations to the unique spin dynamics in geometrically 
frustrated spins. The Curie- Weiss temperature is be- 
tween 61 K to 55 K (using the coupling constants mea- 
sured in the 7 phase), larger than the Neel tempera- 
ture Tn of 35 K. The frustration parameter, defined as 
|0CwH/Ttv = 2.3 ~ 2.5, is large. In temperatures be- 
tween Tjv and |6cw|) th e quantum fluctuation of the 
spin could lead to the spin liquid state. The formation 
of the quantum spin liquid requires dominant direct ex- 
change, that can be described by Heisenbcrg Hamilto- 
nian. The intrinsic fluctuations lead to a small coherence 
scale. In addition to quantum fluctuations, the spin-orbit 
coupling and inter-plane coupling can both give rise to 
fluctuations, lifting the degeneracy of the ground state. 
The bandwidth of the energy excitation is on the order 
of 6cw| = 4.7 ~ 5.3 meV, which agrees with the fit- 
ted values of T (Fig. [7b) and is too big to be accounted 
for by thermal fluctuations alone. As an exercise, we 
could attempt to extract the scaling law of this transi- 
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FIG. 8: The temperature evolution of the magnetic dominant 
S(Q,u>) (integrated over Q = 1.2 ~ 1.3A _1 ). 

tion by fitting \' with a power law along both sides of 
the transition temperature (Fig. [7b). The critical expo- 
nent is about 1/2 when approaching the transition from 
above, and about -1/3 when approaching the transition 
from below. 

The Q dependence of the elastic scattering intensity 
(Fig. [3J2) is best described by a short range, two dimen- 
sional helicoidal order with a correlation length of 5 A, 
corresponding to the next nearest neighbor on the ab 
basal plane. The three-sublattice Yaffet-Kittel model is 
thus ruled out. The dip on the inelastic scattering in- 
tensity shown in Fig. [3^ (summed over 2~9 meV) is an 
artifact of a dead detector, otherwise, it shows the broad 
peak centering around Q = 1.3 A . The Q width of 
the elastic magnetic peak decreases by about 20% as the 
j3-0 2 cools from 44 K to 24 K. 

Due to the complexity of the spin dynamics in /3-02, we 
don't yet have an energy dispersion relation with which 
to construct a complete map of S(Q,uj), like what we 
have done for a-C>2. 



V. EVOLUTION OF THE DYNAMICS 

Up to now, we have limited our discussions on the 
magnetic excitations, which arc dominating the dynamics 
for Q below 1.5A within the kinematic range covered 
by our INS measurements. Even without the polariza- 
tion analysis, we can safely isolate the magnetic compo- 
nents by taking advantage of the similar lattice structure 



and lattice spacings between a and /3-02, which war- 
rant the overall similarity of the translational excitations 
(phonon) between these two phases. To demonstrate 
this, we plot the dispersion curves of acoustic phonon 
branches along major symmetry axes on the S(Q,u>) 
maps of a and /3-Oa (Fig. [3^, b). These phone disper- 
sion curves clearly outline the regions where the phonon 
excitations dominate. In particular, within the second 
Brillouin zone, one can clearly see the typical excitations 
of acoustic phonons extending all the way down to zero 
lj at Q = 2.2 ~ 2.4A , connecting to the elastic nu- 
clear Bragg lines (Fig. [2^, b). The phonon excitations 
within the first Brillouin zone are out of the assessable 
kinematic range. At larger Qs, the multiplicities of the 
inverse lattice planes increases, and it becomes harder to 
identify the specific phonon branches as they overlap. 

Furthermore, the magnon and phonon dynamics can 
be well separated because the molecular spin form factor 
suppresses the scattering intensity at large Q, limiting 
the magnetic excitations to Q < 2A . To study the 
evolution of spin dynamics over solid phases, we integrate 

S(Q,oj) over Q = 1.2 ~ 1.3A , and plot the scattering 
intensity as a function of ui and temperature (Fig. [5]) . As 
the solid cools through the j3-a transition, the magnetic 
excitation transforms from the broad-band dispersion- 
less excitations into the well defined magnon excitations 
consisting of a high energy mode around lOmeV and a 
low energy mode around 6 meV, leaving almost no scat- 
tering amplitudes at low ui < 4 meV in the a phase. The 
same details are illustrated in Fig. [9jt, b comparing the 
scattering amplitudes between a-02 at 4 K and at 
27 K. 

The evolution of the phonon dynamics can be exam- 
ined by plotting S(Q,lo) at constant Q slices at Q = 

1.9 ~ 2.5A (Fig. IHt, d). Solid P-O2 shows more quasi- 
clastic scattering around the phonon branches when com- 
pared to the a phase. This might be a result of the 
additional spin-orbit and intcr-plane couplings in the /3 
phase transferring the high degree of spin fluctuations 
to the lattice vibration. Beyond the additional quasi- 
clastic scattering amplitude, there are not much notice- 
able differences in the phonon dynamics between these 
two phases. 

We may gain insights into how the dynamics evolve by 
looking into the density of states. Dividing the dynamic 
structure function S(Q, ui) by the trivial temperature de- 
pendent factors, including the Dcbyc- Waller factor e~ 2W 
and the phonon (or magnon) occupation number n(u), T), 
and integrating over the appropriate Q range allows us 
to extract the density of states of excitation modes of 
interest. 

fQmax 2W 

J Q . h' ! Q z /2m (n(u),T) + 1) 

(16) 

where [QmimQmax] = [1.0, 1.5]A for magnetic exci- 
tations and [QmimQmax] = [1.7, 5.0]A for phonons. 
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FIG. 9: Constant Q slices of S(Q,lj in a-02 at 4 K and P-O2 
at 27 K. Magnetic excitations (a, b) and phonon excitations 
through nuclear scattering (c, d). 



The resulting density of states are plotted in Fig. [TO] 
The phonon density of states in CY-O2 can be fitted by a 
quadratic function w 2 , as expected, and the magnon den- 
sity of states in CV-O2 follows the theoretical predictions 
as discussed in sec. IIV Al For /3-02, it is interesting to 
note that both nuclear and magnetic density of states can 
be fitted by a linear function in cj with the same slope, 
indicating that the fluctuations arc correlated and might 
be due to the same cause. The density of states for 7-O2 
at high Q is larger than that of /3-O2 as a result of the 
additional orbital rotational diffusion 2 ^. 



VI. UCN CROSS SECTIONS 

The energy of UCN (~ 100 neV) is several orders of 
magnitudes smaller than the energy scale of cold neutrons 
(~ few meV) . For scatterings involving UCN production, 
because of the small momentum carried by UCN, k ucn , 
the momentum transfer Q is approximately equal to the 
initial momentum of the incident neutrons. Similarly, 
Q is equal to the final momentum of the upscattered 
neutrons in UCN upscattering. Therefore, the differential 
cross sections for UCN production (via downscattering) 
and UCN upscattering can be attained by simply reading 
the value of S(Q, u> = Ef ree ) under the parabolic curve of 
free neutron dispersion, Ef ree = h 2 Q 2 /2m, as illustrated 
by the solid red curve (for downscattering) and dashed 
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FIG. 10: The density of states of magnon and phonons in s- 
O2 (a) total excitations (b) magnetic component (c) phonon 
component. 



curve (for upscattering) in Fig. [3jt, b, c, and Fig. [5] 
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At this stage, no effort was made to obtain the absolute 
cross section on the INS data, and the UCN cross sections 
are presented on a relative scale with an arbitrary nor- 
malization. These temperature-dependent cross sections 
are then used as essential input parameters in the Monte- 
Carlo simulations to generate UCN yields that can be 
directly compared to the experimental UCN production 
data—. If desired, the result of the comparison through 
a x 2 analysis yields the overall normalization required to 
set the absolute scale of the cross sections, but we omit 
this step here. 

We plot the resulting (relative) differential cross sec- 
tions for UCN production (Fig. ITTa ) and UCN upscat- 
tering (Fig. [TTb ) in S-O2 in the a phase at 4 K, ft phase 
at 25 K and 42 K, and 7 phase at 47 K. The differential 
cross section for UCN production can be integrated with 
the incident cold neutron spectrum to attain the UCN 
production rate specific to each experiment. Note that 
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FIG. 11: (a) Differential cross sections for UCN production 
in S-O2. (b) UCN upscattering cross sections in S-O2. Note 
that in both plots the low E data (< 1 meV) is contaminated 
by the finite width of elastic scattering. 

the production rate can be maximized by choosing an ap- 
propriate cold neutron spectrum that optimally overlaps 
with the differential cross section — ■ As illustrated 
in Fig. [TTr . the /3-C>2 prefers a colder CN spectrum for 
maximum UCN production. The UCN upscattering, on 
the other hand, is not influenced by the CN spectrum, 
but instead is directly proportional to the number of 
magnons (or phonons) present in the solid and increases 
exponentially with increasing temperatures. To estimate 
the total UCN upscattering rate at different tempera- 
tures, we integrate Eg . IT51 over all final phase space. 

A totally thermalized flux can be described by a 
Maxwcll-Boltzmann distribution characterized by the 
mean temperature of the energy spectrum. With a 
Maxwellian CN spectrum of temperature Tew, we can 
calculate the UCN production following the prescription 
described above to study the dependence of UCN pro- 
duction on the CN spectrum. The resulting temperature- 
dependent UCN production rates for S-O2 coupled to CN 
flux of several Tcn are shown in Fig. [T3] Depending 
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FIG. 12: Temperature-dependent UCN production rates in s- 
O2, coupled to CN flux with a Maxwell-Boltzmann spectrum 
of 10 K, 20 K, 30 K and 60 K. 

on the Tcn, the UCN production in a-0 2 is 2~4 times 
smaller than that in the /3-02- The dominant mechanism 
of UCN production in a-02 is probably through phonon 
coupling around 6 meV, where the free neutron disper- 
sion curve intersects with the phonon branch (as shown in 
Fig. EJi) . Unfortunately, due to the large spin anisotropy 
energy, the kinematical range where the magnon exci- 
tation becomes strong lies above and totally misses the 
free neutron dispersion curve. The magnon scattering, 
even though strong in the a-02, is inaccessible to the 
inelastic scattering process that produces UCN through 
the supcrthcrmal principle. Above the phase transition, 
the spin dynamics changes into the highly unstable state 
with strong short-range correlations as a result of the ge- 
ometrical frustration. The resulting dispersion-less mag- 
netic excitations spread over u, extending down to w = 0. 
These magnetic excitations overlap with the free neu- 
tron dispersion curve from Q = 1.0 ~ 1.7 A , enhancing 
the UCN production for CN with energy ranging from 
2-6 meV (as illustrated in Fig. [Tljt)- The UCN pro- 
duction in 7-O2 is even larger, but the correspondingly 
large upscattering rate makes the UCN loss too big to be 
useful. 



VII. CONCLUSION 

In this work, we update the seminal work of Stephens 
and coworkers^ using the modern high resolution INS 
measurements. Our INS data confirms the detailed pre- 
diction of the 2D AF magnon model, and we report new 
values of anisotropy energies that prefer the 2D Ising 
model over the XY model suggested by Stephens et al. 
Our initial over-estimate of the UCN production in a-02 
is a result of the low values of spin anisotropy energy sug- 
gested by Stephens et al. With the corrected anisotropy 
energies, we confirm that the UCN production in a-02 is 
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predominantly through phonon coupling, and thus suffers 
from the mass suppression factor because it is harder for a 
neutron to nuclear recoil the whole O2 molecule that ini- 
tiates the translational wave. In the /3 phase, we have re- 
ported an unexpected large UCN production^. The INS 
measurements unfold the origin of the enhanced UCN 
production in /3-02 as the unique dispersion-less mag- 
netic excitations observed in other similar geometrically 
frustrated spin systems ^ 10 ' 25 . The spin dynamic in (3-0 2 
has been previously characterized using a phenomcnolog- 
ical model of coupled paramagnets^ on a three sublat- 
tice Yaffet-Kittel structure^, however, the ground state 



configuration must be highly unstable and the INS data 
shows many distinct features for the intriguing possibil- 
ity for spinon excitations. The description of the detailed 
dynamics still awaits for better INS data using single 
crystal samples, which are challenging to prepare due to 
the thermal properties of oxygen. 
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